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Abstract 

Experimental data for bb quarkonia have been compared with the 
predictions of a variety of nonrelativistic quark models. It is found 
that a potential a^Jr — b/r+ const gives good agreement, while many 
others do not. Some implications of this observation are discussed. 



1 Introduction 



According to a recent review the masses of the hh quarkonia below twice 
5^ \ the mass of the B meson are successfully predicted by a variety of potential 

models. The average error of the predictions is typically a few MeV — 2.3 
MeV in the best case[0. Some of these models include relativistic correc- 
tions, but if we replace the masses of the x multiplets by their centres of 
gravity, the relativistic corrections are not crucial for the success. This re- 
sult is puzzling. From fine splittings, which are a purely relativistic effect, 
or from the differences between typical kinetic energies evaluated relativisti- 
cally and nonrelativistically, one finds that the relativistic corrections should 
be of some tens of MeV. The standard explanation is that what is called 
the nonrelativistic Hamiltonian, is in fact an effective Hamiltonian with the 
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parameters modified by relativistic corrections. This rises the question: is it 
possible to reproduce the data within errors using such an effective Hamilto- 
nian? Since experimental errors on the masses are of the order of 0.3 MeV, 
even the " successful" models correspond to very low confidence levels as eval- 
uated from the test. Thus in order to answer the question one should do 
better. Moreover, there are also other data, e.g. leptonic widths and electric 
dipole transition probabilities, which could and should be fitted. Quark - 
antiquark potentials derived from studies of quarkonia are used in a vari- 
ety of applications. Let us mention as examples applications to heavy light 
systems^, to be mesons[^, and to ti production^. Therefore, improving on 
existing potentials seems useful. On the other hand, as is well known, two 
body relativistic quantum mechanics does not exist yet. If no nonrelativistic 
potential can reproduce the data, this may be a valuable hint on how to 
construct a relativistic theory. 



2 The model 

There is no standard nonrelativistic quarkonium model. The Schrodinger 
equation is 

- —V^iP + V(rU = EiP, (1) 
rrib 

but the quark mass rUb and the potential V{r) vary from paper to paper. 
Fortunately, most proposals for the potential either are of the form 

V{r) = -ar-"" + br'^ + Ct, (2) 

where a,b,a, (3,Ct are nonnegative constants, or are numerically very well 
approximated by such formulae. Let us mention as examples the Cornell 
potential^ with a = (3 = 1, the potential of Lichtenberg and collaborators 
with a = P = 0.75, the potential of Song and Liu^] with a = (3 = 0.5, the 
logarithmic potential of Quigg and Rosner^ corresponding to a = 
and the potential of Martin|[T^ corresponding to a = 0, (3 = 0.1. Other 



potentials, which have sometimes quite different analytic forms, become very 
close to potentials of this type, when their parameters are adjusted to fit 
the data. We have checked that in particular for the very successful Indiana 
potential [jlT| and for the famous Richardson potential|T^. We found that 
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it is enough to impose on all these potentials the condition that they fit 
the energy spectrum as well as possible within the freedom described below. 
This is related to the suggestion of Eichten and Quigg, who concluded from 
their inverse scattering analysis [|13|] that the first few L = energy levels 
and the corresponding leptonic widths (not used by us, but we use instead 
the energies of the L = 1 levels) to a large extent determine the potential in 
the region relevant for the quarkonium calculations. Note, however, that we 
have no proof that a completely different potential would not fit the data as 
well. 

In order to make our search for a satisfactory potential more effective, we 
introduced a scaled variable p = Xr. In terms of this variable the Schrodinger 
equation (1) with potential (2) can be reduced to the form 

-V^^(p) + V(p)^(p)=^^(p), (3) 

where 




(4) 
(5) 

"^(E-Ct), (6) 

(7) 

Thus, for given values of the parameters a, (3 the model has four free param- 
eters, which can be chosen as ruh, A, Ct, and C. Our idea is to choose a set of 
observables, which depend only on the parameter C. Parameter C is chosen 
so as to minimize fo'^ the observables chosen. Then other observables can 
be used to determine the remaining parameters of the model, but this does 
not influence the quality of the fit. Using this procedure for the potentials 
proposed in the literature we got significant improvements of the fits, but 
none was quite satisfactory. 



3 Observables 

As observables we have chosen 
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6, = 0.6290 ±0.0005. (S) 



M(35) - M{1S) 
M{3S) - M{2P) 
M{2S) - M(1P) 
M(2S) - M(1P) 
M(25) - M(15) 

IV'25(0)P 



0.774 ±0.006, (9) 
0.219 ±0.001, (10) 



0.492 ±0.111, (11) 
0.433 ±0.071, (12) 



I^15(0)P 

k = |V^is(0)|i(lP|r|25) = 2.29 ±0.16, (13) 
= |^/;i5(0)|i(2P|r|35) = 1.59 ±0.15, (14) 

The observables bi,...,bs depend on the parameter C, but not on the 
other three parameters. Thus the distribution corresponds to seven de- 
grees of freedom. All the numerical values are calculated from the data given 



in the 1994 Particle Data Group Tables [|I4|. There are some problems with 
the observables feg and 67, however. The formulae connecting the measured 
leptonic widths to the calculated |'?/'ri,s'(0)P contain a radiative correction fac- 
tor. This is 1 — and we put it equal 0.7. Since neither the exact value of 
as, nor the effect of the higher order terms is known, this number is uncertain 
by perhaps some 20%.. This uncertainty very probably cancels in 64, 65 and 
bs , but for bg and 67 it introduces a common factor, which deviates from 
unity by some seven per cent. This uncertainty has not been included in the 
errors. Note also that the simple dipole formulae are much less reliable for 
quarkonia than for atoms. 

Once the parameter C is chosen so as to reproduce as well as possible 
the observables &i, . . . ,68, the parameters A, mb and Ct can be calculated 
from the experimental value the leptonic width of the T(IS') state, which 
yields |?/'is'(0)p, from the mass difference M{3S) — M{1S) and from the 
mass M[1S). This calculation has no bearing on the quality of the model. 
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except that perhaps values of mf, too far from 5 GeV would make the model 
unplausible. 



4 Results and conclusions 

After some exploration in the a, (3 plane, we have found a positive answer to 

our question. For « = 1, (3 = 0.5 we find for C = 0.9315 the value = 6.5, 
which is very satisfactory for the seven degrees of freedom. This solution 
corresponds to the potential in eq. (1) given by 



where V{r) and r are in GeV. The corresponding quark mass mf, = 4.79333 
GeV is quite reasonable. For r — > our potential has the r^^ dependence 
corresponding to one gluon exchange. With present data, however, we have 
no evidence for the additional factor l/log(Ar), which according to QCD 
should be introduced by the running of the coupling constant. The expected 
part of the potential linear in r is not seen. Probably the bottomonia are too 
small to reach sufficiently far into the asymptotic region of linear confinement. 
Perhaps a more flexible potential would exhibit the linear part. 
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